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ABSTRACT 



Discontinuous initial wave functions or wave functions with discontintuous derivative 
qh and with bounded support arise in a natural way in various situations in physics, in 

particular in measurement theory. The propagation of such initial wave functions is not 
well described by the Schrodinger current which vanishes on the boundary of the support 
of the wave function. This propagation gives rise to a uni-directional current at the 
boundary of the support. We use path integrals to define current and uni-directional 
current and give a direct derivation of the expression for current from the path integral 
formulation for both diffusion and quantum mechanics. Furthermore, we give an explicit 
asymptotic expression for the short time propagation of initial wave function with compact 
support for both the cases of discontinuous derivative and discontinuous wave function. 
We show that in the former case the probability propagated across the boundary of the 
support in time At is O (At 3 / 2 ) and the initial uni-directional current is O [At 1 / 2 ). This 
recovers the Zeno effect for continuous detection of a particle in a given domain. For 
the latter case the probability propagated across the boundary of the support in time At 
is O (At 1 / 2 ) and the initial uni-directional current is O (At- 1 / 2 ). This is an anti-Zeno 
effect. However, the probability propagated across a point located at a finite distance 
from the boundary of the support is O {At) . This gives a decay law. 



1. Introduction 

Detection of the coordinate of a quantum particle can give the negative result that 
the particle is not in a given domain, for example, it is not in a half space bounded by 
a given plane. This means that at the moment of measurement the wave function of 
the undetected particle vanishes identically in the given domain, for example, beyond the 
separating plane. It follows that the wave function of the particle is supported outside 
the measured domain. The boundary of the support is the boundary of the domain, for 
example, it can be the plane separating the measured half space from the other half. The 
instantaneous killing of the wave function in a domain G for a given time interval [to, to + 
At] can be accomplished mathematically by introducing in the Schrodinger equation the 
(dimensionless) time dependent potential 

V(x,t) = OO • XG(x)x[t ,t +At](t) 

where 

f 1 if x G G 
Xg{x) - I q otherwise> 

After the instantaneous negative detection the particle in the time interval [to, to + At] it 
is free to propagate into the domain, that is, the potential V(x,t) is turned off instanta- 
neously. A discontinuous wave function with compact support can be realized mathemat- 
ically by introducing a similar potential that contains the spatial derivative of the delta 
function on te boundary of the domain and an infinite potential inside the domain. 

The propagation of the wave function of the particle after the negative measurement 
is not well described by the Schrodinger current because the compactly supported initial 
wave function after the measurement gives rise to an initially vanishing Schrodinger cur- 
rent at the boundary of the support. Indeed, assuming the wave function is continuous, 
since 

ijj(x, t)\ dG = 0, for all t < t < t + At, 

where t is the instant of the measurement and dG is the boundary of the measured 
domain, we must have 

J dG (x,t) = —Qmip(x,tWip(x,t) = 0, for all t < t < t + At, 
m 

It is obvious, however, that despite the vanishing current on dG there is propagation 
across dG after the instant t + At. This is in general the case in Schrodinger's equation 
if the initial wave function has compact support. The case of a discontinuous initial wave 
function is problematic as well. 
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It is the purpose of this paper to calculate the short time propagation of the wave 
function across the boundary of the support of the initial wave function. This propa- 
gation gives rise to an instantaneous uni-directional current into the measured domain 
G. Our main result is an explicit asymptotic expression for the wave function at time 
to + At for small At. Our analysis is one- dimensional, however the generalization to higher 
dimensions is straightforward. 

/,From this expression, we find that for a continuous wave function the probability that 
propagates across dG in time At is O (At 3 / 2 ) and the rate of propagation is O (At 1 / 2 ). 

The latter means that the initial uni-directional current is O (At 1 / 2 ). We also determine 
the probability P c that propagates in time At beyond a distance c from dG, 



P 



c 




This expansion is valid for c > yfa. For example, if c = O (At 1 / 3 ') , we obtain P c = 

0(At 2 ). For c = 0(VAt) the propagated probability is 0(At 3 / 2 ). This gives an 
estimate on the dependence of the probability of detection on the size of the detector. 
Detection at a point or in a domain is discussed in the context of absorption in separate 
papers gj-[f§. 

Our result recovers the Zeno effect [||, |5| that a quantum particle cannot arrive at a 
point under continuous observation (the particle is "frozen" or "reflected back" by the 
continuous observation). This result is in agreement with the result of [0] (see Section 6 
for more detailed discussion). 

For a discontinuous initial wave function, we find that the probability that propagates 
across dG in time At is O (At 1 / 2 ) and the rate of propagation is O (At x / 2 ) . The latter 

means that the initial uni-directional current is O (At- 1 ' 2 ) and is initially infinite. This 
result is an anti-Zeno effect M. The probability P c for this case is given by 

This implies that the survival probability decays exponentially (see eq. (|5.1| )). 

Wave functions with compact support and discontinuous derivatives appear in other 
applications as well 0, pj. They represent actual physical situations and are of both 
theoretical and practical interest. 

The notion of uni-directional current can be used to define the notion of "time of 
arrival" and measurement in quantum mechanics. A review of different approaches to 
the definition of uni-directional current related to time of arrival in quantum mechanics 



is given in |10|, |TT| . Our approach to the definition of a uni-directional current is different 
from all the above mentioned approaches. It is defined only for wave functions that vanish 
identically beyond the detector, that is, its support is initially bounded by the detector 
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and is denned only at the boundary of the support. At such points the uni-directional 
current and the current are the same. 

Out approach to the calculation of the uni-directional current is based on a Feynman 
integral approach. This seems to be a new approach that is fundamentally different 
than the above mentioned approaches in that it is not approximate and gives the full 
asymptotic behavior of the wave function for short times. A direct derivation of the 
Schrodinger current from the Feynman integral is given and extended to the case of a 
uni-directional current. 

In this approach uni-directional currents appear when the wave function satisfies the 
following conditions, 

• The initial wave function has compact support. 

• At the boundary of the support either the wave function or its derivative has a 
discontinuity. 

Initial wave functions that satisfy these conditions are denoted ICSWF (initial compact 
support wave function). In general, the solutions of Schrodinger's equation do not develop 
discontinuities or non-smoothness. However, there are situations where non-smoothness 
does occur 0. If an infinite potential barrier is introduced, the gradient of the wave 
function has a discontinuity across the barrier. Indeed, in our formalism for the description 
of an absorbing wall such a discontinuity arises in a natural way. 

In Sections 2, the case of compact support in diffusion theory is reviewed. The Wiener 
path integral formulation of diffusion is shown to lead to the classical definitions of dif- 
fusion current and to a uni-directional diffusion current at an absorbing boundary. This 
derivation seems to be new. In Section 3, a Feynman integral formulation is used to define 
quantum probability current and a new direct derivation of the Schrodinger current from 
Feynman's integral is given. This definition leads to the definition of a uni-directional 
current at the boundary of the support of the initial wave function. In Section 4, the 
main results of the paper are presented. These are short time asymptotic expansions for 
the propagation of an ICSWF and the uni-directional current. Finally, Section 6 contains 
a discussion and summary of the results. 

2. Current and uni-directional current in diffusion 



We consider a diffusion process, x(t), with noise coefficient a(x,t) and drift b(x,t) 
|T2|| . The transition probability density of the process is denoted p(y, t\x,s). This is the 



probability density 

d 

p(y, 1 1 x, s) = — Pr {x(t) < y \ x(s) = x] 
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It satisfies the Fokker-Planck equation 







1 d 2 



m P(y,t\x,s) 2dy2 



d_ 

dy 



with the initial condition 



limp(y, t 1 x, s) — 5{y — x). 

tls 



(2.1) 



(2.2) 



At an absorbing boundary at y = 0, say, the function p(y,t \ x,s) satisfies the boundary 
condition 



p(0, t 1 x, s) = 0. 



(2.3) 



The uni-directional probability current (flux) density at a point X\ is the probability 
per unit time of diffusing trajectories that propagate from the ray x < X\ into the ray 
x > X\. It is therefore given by 



Jlr{xi, t) = \im Q J L r{xi, t, At), 



where 



JLR(Xl,t,At) 



(2.4) 



(2.5) 



dy 



_ / dx / — 

At Jxi J-oo v / 27rAt a(y, t) 



<J X - V -Hy,t)*tf} p( M 

'1 2a ! (t,,()A( J" 9 ' 



(note that the dependence of p on the backward variables x, s has been suppressed). The 
integral (p75| ) can be calculated by the Laplace method at the saddle point x — y = X\. 
First, we change variables in (|2.5| ) to 



to obtain 



JLR(Xl,t,At) 



x = X\ + £ v At, y = x\ — f]v At 

dr\ 



o V2nAta(x 1 - r]VAt,t) 



exp < 



£ + 7] - b{x x - 77V At, t) V At 



2a 2 (xi - r^VAt, t) 



p (xi — r/VAt, t — At 



and changing the variable in the inner integral to 77 = ( — £, we get 

-(C-0 Vam) 



/■OO /"OO 

J LR ( Xl ,t,At)= dU r—- , 

Jo J£ v27rAtcr(xi 



x 



(2.6) 



exp 



c-^i-(c-e) vAt,t)VAt 



2a 2 (x 1 -(C-0 VAt,t) 



pfci - (C - v 7 ^, t - At 



Next, we expand the exponent in powers of v At to obtain 



2a 2 (x!-(C-0 VAM) 

C 2 , \<?{C-t)o*(xut) CK^nt) 



(2.7) 



2<x 2 (xM) 
the pre-exponential factor, 
1 



2a 4 (x 1 ,t) 



At + O (At) 



<t (a* - (C - V AM 
and the pdf 



a(xi,t) 



l + ^ 1 4(C-0v / At + O(At) 
a(x ± ,t) 



(2.8) 



p(^i-(C-0 VAt,t-At 
P(x 1 ,t)- 9p ^ t) (C-QVAt 



(2.9) 



Using the expansions 

J LR (x 1 ,t,At) 



OO /'OO 



<9x <9t 
79T) in (12.9), we obtain 



£At + O (At 3/2 ) 



(2.10) 



■ exp 



dp(x 1 ,t) 



-VAt 



{ v / 2S<r(xi,t) 1 2a 2 (a;i,t)J f # 
YC 2 (C - J {xi,t) {b(x u t) a' ( Xl ,t) (C - 



At 



cr 3 (xi,t) 



a 2 (x h t) 



a(x±,t) 



p(xi,t) 



+ (C-£)p' (*i,t)]+0(At)}. 



Similarly, 
where 



J RL (x l ,t, At) 
1 z -21 



JnL(xi,t) = Km JRL{xi,t, At), 



,/,, ^ ^ Jx-y-b(y,t)Atn ^ t _ My 



(2.11) 



2a 2 (y,t)At 



At, — OO JI1 V27rA7o-(?/,i) 
The change of variables in Q2.ll ) 

x = X\ — ^y/At, y = xi + rjVAt 



5 



nves 



JjiL(xi,t,At) 



(2.12) 



■ exp , 

n J* ^/2itAta(x l ,t) [ 2a 2 (x 1 ,t) # 

/C 2 (C - r'Orx, *) C6(^i. *) (C - 



c 2 



+VAt 



a 3 (xi,t) 

(c-OpW) 



cr 2 (xi,t) 



a(xi,t) 



dt 
p(xi,t) 



+ 0(At)\. 



Since p(xi,t) > 0, both Jlr{x±, t) and Jrl{xi, t) are infinite, however, the net flux density 
is finite and is given by 

J net( x i^) = |im o {J LR (xi,t, At) - J RL {xi, t, At)} = 



poo roc fa 

-2 / dU 

A) Jt v / 2vrAta(xi,t) 



exp 



2cr 2 ( Xll t) J 

C 2 (C - ^(zi, *) C&Oi, *) ^i, *) (C - ' 



<J 3 (xi,t) cr 2 (xi,t) 
9 a 2 (x,t) 
~ ~dx~ 2 



cr(xi,t) 
t) + b(x, t)p(x, t) 



P(xi,t) + (C-0P'(^1,<) 

(2.13) 



In deriving eq. ( |2.13|) use has been made of the following identities, obtained by changing 
the order of integration, 



Si j£ 



d£ / 



27T <J 



d£ / 

J£ 



(C - o 



2ti a 



exp 



exp 



2a 2 
2a 2 

"2 ' 



C 4 ^C 

(i 2v / 2tt(T 

- ( 2 d( 



exp 



c 2 



3a 4 



2a 2 



c 1 



2a 2 



V27T(T 

1 C^C 

2 io 



27T <J 



exp 



2(7 2 



(7 

T 



exp 



2a 2 



a 

T 



Equation ( |2.13| ) is the classical expression for the probability (or heat) current in diffusion 
theory [Q. The Fokker-Planck equation (|2.1j) can be written in terms of the flux density 
function J(x, t) in the conservation law form 



§- t p(x,t) = -^-J(x,t). 



(2.14) 



Next, we calculate the uni-directional flux at the absorbing boundary x — 0. The 
absorbing boundary condition (|2.3| ) implies that the pdf vanishes for all x > so that its 
right derivatives at the origin vanish. It follows from eq.( |2.11 ) that 



J RL {o,t) = o. 
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On the other hand, eqs. (|2.5| ) and ( |2.13| ) give 



J(0,t) = J LR (0,t) 



d a 2 (x,t) 
dx 2 



p{x,t) 



x=0 



Since a 2 (x, t) > and p(x,t) > for x < 0, but p(0, t) = 0, it Since a 2 (x,t) > and 
p(x,t) > for x < 0, but p(0,t) = 0, it follows that J(0,t) > 0. This means that 
there is positive flux into the absorbing boundary so that the probability of trajectories 
that survive in the region to the left of the absorbing boundary, f®^ p{x, t) dx, must be 
a decreasing function of time. This can be seen directly from eq. Q2.14p by integrating it 
with respect to x over the ray (— oo, 0) and using the fact that lim :r __ > _ 00 J(x, t) = 0, 



d r° 

— / p(x,t)dx = -J(0,t) < 0. 



3. Current and uni-directional current in Feynman integrals 

To keep the calculations simple, we consider a particle with an infinite potential for 
x > and its free propagation after the infinite potential is turned off at time t = 0. There 
is no analog to uni-directional current in Feynman integrals due to the non-additivity of 
probability on sets of trajectories. First, we examine the notion of current in the usual 
Feynman integral. The current at a point is the net rate of change of probability on one 
side of the point. That is, the current at x — 0, say, is 



lim — , 

A^O At J-oo 



J(0,t) 

According to the Feynman formalism, we write 

if)(x,t + At) 

so that 



\i/>(x,t + At)\ 2 - |V0M)l' 



dx. 



(3.1) 



m f 00 , . . ( im 



dy 



[° \if>(x,t+ At) | 2 dx 



m 



2nihAt J- 



ij}(y,t)exp 



( im 
\2hAt 



[x - yf \ dy 



dx. 



Thus the current is the difference between the population of trajectories that propagates 
in the time interval (t,t + At) from the entire line into the interval (— oo,0) and the 
population there at time t. 

Expanding 



ip(y,t) = ip(x,t) + (y- x)ip x (x,t) + \(y - x) 2 ip xx (x, t) + ■ 



(3.2) 
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we obtain 



L III V dL^' t] exp { ^k {x ~ y)2 } dy 

(, , , .,,2 - , , t , .JftAt 



fix = (3.3) 



x,t)\ + ^ei/;(x,t)ip xx (x,t) dx + o(At). (3.4) 

m J 

It follows from integration by parts that 

h r° - h 

J(0) = -3m / ^(z, tU x Jx, t) dx = — Sm^(0, *M X (0, t), (3.5) 
m J-oo m 

because ^(0, ^^(O, i) is real valued. Equation (|3.5| ) is identical to the Schrodinger 
current. This derivation seems to be new. 

Next, we consider the case that the wave function at time t vanishes outside the 
interval [— a, 0]. As mentioned above, this situation arises, for example, if up to time t 
there is an infinite potential for x > and x < —a and the potential for x > is turned 
off instantaneously at time t. This is a mathematical idealization of various physical 
situations (see, e.g., Chapter 6). Since i[>(0,t) = 0, the Schrodinger current at the point 
x = vanishes at time t, according to eq.( |3.5|) . Yet, there is probability flux across x = 0. 
To see this, we evaluate the rate of population change, (|3.1J ). In the case at hand the 
second term in the integrand of ( p.l|) vanishes at time t. We assume first that x = —a 
is a reflecting wall, that is, ip(x,t) = for x < —a for some positive a. Thus there is no 
propagation across x = —a. The case a = oo is considered in Section 5 below. 

The probability of propagating from a given interval [—a, 0] into the ray [0, oo) in the 
time interval (t, t + At), starting with the wave function ip(x, t) in the interval [—a, 0] and 
outside, is given by 

^out= rmy,t + At)\ 2 dy. (3.6) 
Jo 

The wave function at time t + At is given by the free propagator 

The uni-directional current from [—a, 0] into the ray [0, oo) is the current 

J LR (0,t) = Urn — J o \ij(x,t + At)\ 2 dx. (3.8) 

Note that if trajectories are not truncated for x > 0, the wave function at time t + At no 
longer vanishes at x = and so does the Schrodinger current. 
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4. Short time propagation 



To estimate the uni-directional current Jlr{®, t) for short times, we first estimate the 
integral in eq. (|3.8| ) for small At. We begin with an initial wave function i/j(x, 0) that is a 
polynomial 



A? 



Q( x ) = H ( ij 



X' 



3=1 



in the interval [— a, 0], such that Q(— a) = Q(0) = and i/j(x, 0) = otherwise, the free 
propagation from the interval [—a, 0] is given by 

im(x — y) 2 ' 



1)(y,At) 



m 



2irihAt 



Q(x) exp 



2hAt 



dx. 



The boundary condition at the left end of the support of ip(x, 0) is written explicitly as 



N 



3=1 



Setting 



a 



hAt 



m 



the probability mass propagated freely into the positive axis in time At is given by 



\ijj(y,t + At)\ z dy 



Q(x)e i{x - y) /2a dx 



2ira Jo 

We change variables by setting x = y/a£, y = \/ar], £ = ( + r\ to get 



dy. 



\^{y,t + At)\ 2 dy 



1/2 



2tx Jo 

First, we evaluate the inner integral, 



-v 

a/y/a—rj 



Q (v / "(C + V)) d( 



drj. (4.1) 



N 



j n(v) = E<?jVc^ 

3=1 



a/^/a—r] 



(( + V ye^/ 2 d( 



Integration by parts gives 
In(v) = 



N 

3=1 
N 

3=1 
N 

i=i 



-v ^de< 2 / 2 

ir _ (C + vY—— 



- a \ j e i ( a /v / S+'7) /2 

y/a J a/y/a + r] 



N 



3=1 

e ,c 2 /2 ( iiC + vf^ _ (C + vY 

a/y/a-V V C C 2 



(4.2 

JCV2 d (C + vY 



-a/y/a—r) 



d( C 



d( 
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because 



For j — 1, we obtain 



The function 



N 



5>v^(_|) = 5>(-a)' = o. 



3=1 



-iqwfar] 



3=1 



»7 



a/y/a-r) C 2 



is square integrable and its integral is independent of a to leading order. Indeed, the 
integral 

•A2 /n 
-JJ gl 



exists and near 77 = it is bounded by 

-r, g 



c 2 



< 



-vd( = 1 

- C 2 17 



so that 



g^/2 



-a/y/a-f] C 2 

For large 77, we have the asymptotic limit 



dc 



< 77 



dC 



1. 



-r, e < 2 /2 _ gV/2 



-d( ~ i- 



3 C 77° 

as is easily seen from l'Hospitale's rule, so that 



n 



e*C 2 /2 



%/y/a-T] C 2 



dC 



1 

< — . 

- ^4 



This means that the function 



V S-a/y/a-n 6 C 2 ^ 



bution to the integral ( |4.1|) is to leading order 

a 3/2 , 



is integrable. It follows that its contri- 



2tt 



I9i I 



V 



v e i C 2 /2 2 



c 2 



dC 



(4.3) 



Now, we consider the term j = 2 : 



iq 2 a 



-a/y/a—r) 
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Proceeding as above, we find that both terms in the integral are uniformly square inte- 
grable functions of 77 for all a > sufficiently small. It follows that the contribution of the 
term j = 2 to the integral (|4.1|) is O (a 5 / 2 J. The mixed term involving qiq 2 contributes 
0(a 2 ). 

Next, we consider the third order term: 



-a/y/a-rj C, 



a 



2 " - , 2 2(-V 

(C + v) 



a \ * 2 (o/y/a + r])-r] j( a/ ^ +v ) 2 /2 [-1 j?/2± 
(—a/y/a — rj) J-a/^- v dQ 



O(-). (4.4) 



This term has a pre-factor of q^a 3 / 2 so that its contribution to ( |4.1| ) is O ^a 3 / 2 ^) . Proceed- 
ing by induction, we find that all terms contribute O (a 3//2 ) to ([4. 1|) . It follows that 



h \^(y,t + At)\ 2 dy = 0\^-^j |. (4.5) 

If a = 00 and all Fresnel-type integrals are interpreted as the limits 

f(x)e ix2/2 dx = lim [ V f(x)e { - e+i)x2/2 dx, 



fhAt\ 



3/2N 



00 



the term O (1/a) in eq.( O ) is replaced by 0(1). It follows that the higher order terms 
of the polynomial contribute higher order terms in the expansion of In in powers of a. 

Obviously, if the polynomial is replaced by an analytic function that vanishes at the 
ends of the interval, the result remains unchanged. Furthermore, if Q(x) is a square 
integrable analytic function on the negative axis, eq. (|4.5| ) holds. 

The asymptotic estimate (|4.5|) is valid when higher order terms can be neglected 
relative to lower order terms. To get an explicit bound on At from this condition, we 
write 4>(x,t) as a series of eigenf unctions 

iE n t 

ipj{x) exp < - 



^,0 = E^(^)exp{-^}, (4.6) 



where ip n (x) are eigenfunctions that satisfy the boundary condition ^> n (0) = and E n are 
the corresponding eigen energies, to find that 

dxi 3W) 
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and 







—2mE n 


3 






Q x 2j+1 




h 2 




dx 




<9 2j >„(0) 




-2mE n 


3 


^(0)| 




dx 2 i 




h 2 





(4.7) 



of qj in the expansion is O y (j^j 
the expansion is that 



The asymptotic evaluation of the integrals for ?/ near or at gives that the coefficient 

. It follows that the condition for the validity of 

Q2j+l(t) 



hAt 



III 



« ^7t( (4-8) 

for all j > 0. Using eqs.(4.6)- (|4.7Q in ( |4.8|) , we obtain that the condition for the validity 
of the expansion is 



hAt 



m 



(2i + l)(2j + 3)E^i(=^) 



2mE n y dipnjO) -jEnt/h 



dx 



v «d ( ~2mE n \ j+1 dip n (0) iE n t/h 
^n=l { h 2 ) dx e 



(4.9) 



for all j > 1. 

If, for example, the initial wave function is a single eigenfunction, the condition (|4.9|) 
reduces to 



E n At < h. 



(4.10) 



The analysis of the continuous spectrum case is identical. The summation with respect 
to n in the condition (|4.9Q is replaced by integration with respect to n. 

It follows from eq. (|4.5| ) that for short times 

JlrM = 0(Vt) 

so that the population in y > increases as 0(t 3 ^ 2 ) for short times. Obviously, once 
if)(Q,t) 7^ 0, the current becomes the usual Schrodinger current. 

The probability mass that propagates in time At beyond a fixed point c > can be 
found from the above expansions. This probability is defined as 

2 



\4>(y,t + At)\ 2 dy 



a 



1/2 



1lX Jc/^/a 



a/^fa—ri 



Q (v / "(C + i?)) e* 2/2 dC, 



d-q, (4.ii) 



rather than (|4. 1|) . The individual terms in the expansion 
with the obvious changes. Inequality ( |4.3| ) becomes 



a 



3/2 



27T 



i 12 r 


r—r] g*C 2 /2 


m\ / 


v L c 2 dC 


Jc/y/a 



drj = O (a 



are estimated as above 



(4.12) 
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The same estimate applies to all terms in the expansion. It follows that 




(4.13) 



This expansion is valid for c > \fa. For example, if c = O (a 1 / 3 ) , we obtain P c = O 
For c = O (y/a) the result ( |4.5|) is recovered. 

Now, we consider an initially discontinuous wave function 



of 




for x < 
for x > 



where Q(x) = YljLoQj^ vanishes at x — —a and Q (0) ^ 0. First, we consider the 
propagation of the term q . Its contribution to the propagated wave function is 



a 



1/2 



2tt Jo 



Qo 



a/y/a—r] 



d( 



d v = (a 1 ' 2 ) . 



This gives rise to an infinite uni-directional current at the point of discontinuity. Next, 
we calculate the probability propagated beyond a distance c > from the discontinuity. 
Again, the term q contributes 



a 



1/2 



2lT Jc/^/a 



Qo 



i/y/a—rj 



drj 



M 2 « +Q (a z/T 



2nc 



c 



(4.14) 



It follows that for c = O {a 1 / 2 e ) the propagated probability is O (a l / 2+£ ^j , so that the 
resulting uni-directional current is infinite for < e < 1/2. If c = O (1), the propagated 
probability is O (At) and gives rise to a finite current. 

The analysis of the propagation of a polynomial shows that the other terms in the 
polynomial contribute higher order terms to the propagated probability. 

5. Discussion and summary 



The expression ( |4.5|) can be applied to the following experiment. A particle with 
energy E n is released between two perfectly reflecting walls placed at x = and x = —a 
(a > 0). A perfectly absorbing detector is placed at xo > 0. The reflecting wall at x = 
is removed instantaneously for a time interval of length At and is then instantaneously 
reinstated. If the particle crosses in this time interval, it gets registered by the detector. 
According to eq.( |4.5| ), if this experiment is repeated iV times, the number of particles 
registered by the detector will be proportional to iV (At) 3 ^ 2 , if At satisfies the condition 
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Note that the condition ( [4.1 0| ) in this context has no relation to the energy-time 



uncertainty principle because At is not related to a measurement of the particle's energy. 
The removal and reinstatement of the reflecting wall changes the energy of the particle 
so that it is not known what it is after the wall is reinstated. 

A possible physical approximate realization of this experiment consists in turning off 
and on again a detecting device, e.g., by illuminating the detection region (and using e.g., 
the Compton effect), a time interval of length At. According to eq. (|4.5| ), the survival 
probability on the left of the barrier is 



S (At) = 1 - c (At) 



3/2 



where c is a constant. If this experiment is run on N identical systems with NAt = T 
and At satisfies the condition (|4.10 ), the probability that all survive by time At, denoted 



S T (N) = (l-c(At) 

r 

= 1 - 



3/2 



, N 




N 



'N 

The expected number of systems that decay is 

(N) = O (iVAt 3/2 ) = O (TVAt) ^ for N > 1. 



The result ( |4.5|) is similar to that obtained in time dependent perturbation theory, 
known as the Zeno effect H [7j, where the probability that an irreversible decay will 
occur before a short time t is O (t 2 ). As in the Zeno effect, the law ( f4.5| ) indicates that 
propagation into a detection region under continuous observation makes it impossible for 
a particle to be observed. This phenomenon is referred to as the freezing of a particle in 
its initial state. This apparent paradox disappears in quantum theory with a measuring 
device, as shown in 0, |3J]. 

We return now to the ideal detection experiment by illuminating the detection region 
x > Xq for a short time At at time intervals At apart. The detection probability is 
P c , given in eq.( |4.11] ). If x = 0, the detection probability is O (a 3 / 2 ) per illumination 



pulse. The result remains the same if the illuminated region is the interval (0, c) , where 
c = O {^y<y)■ If, however, one illumination pulse covers the region x > and the following 
one covers x > c, where c = 0(1) for a satisfying eq.([4.8p or, equivalently, Q4.1U| ), the 



conditional probability of detecting the particle in the second pulse, given that it was 
not detected in the first one, is O (a 3 ). The former result means that the width of the 
illuminated region has to be at least O {\^a) to achieve the maximal order of magnitude 
of the probability of detection per pulse. 

The result eq. ( |4. 14|) can be applied to the following ideal measurement experiment. 



A discontinuous wave function is created by introducing a potential V(x) = 5'(x) M. If 
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the illumination region is x > 0, the probability of detection per pulse is O {^/<y). If the 
the regions x > and x > c are illuminated alternatively, the conditional probability 
of detecting the particle in the second pulse, given that it was not observed during the 
first one, is 0(a), according to eq. (|4.14j ). This means, that if in the second pulse the 
illumination of the region x > c is kept forever, the conditional survival probability S(t) 
(the probability of not observing the particle by time t after the beginning of the second 
pulse) is 

S(t) = O (expj^jjgo^)! 2 eft'}) (5.1) 

for some 7 > 0. The probability propagated across the boundary of the support in time 
At is O (At 1 / 2 ) and the initial uni-directional current is O (At- 1 / 2 ). This is an anti-Zeno 
effect ||. It means that in N identical systems observed for time At = T/N the the 
probability that all survive by time At is 

S T (N) = (l - c (At) 172 )^ « exp {-cVN} ^ for iV > 1. 
The expected number of systems that decay is 

(N) = O (iVAt 1/2 ) = O (VTN) ^00 for N > 1. 

It should be remarked that the leading order short time asymptotics is unaffected by 
the presence of a finite potential beyond the boundary of the support of the initial wave 
function. This suggests the possibility that the detection region beyond the support of 
the initial wave function can be characterized by a potential without essentially changing 
the above result. 

In summary, we compared the notions of net and uni-directional fluxes in the Wiener 
and Feynman integrals. At points where the density does not vanish the uni-directional 
fluxes are infinite, though the net flux is finite and is given by the traditional expressions 
for flux in the diffusion and Schrodinger equations. At points where the density vanishes, 
for example at certain types of boundaries (absorbing for Wiener trajectories and at the 
boundary of the support of the wave function) the uni-directional fluxes are finite. In 
the Wiener integral the uni-directional flux at an absorbing boundary does not vanish, 
resulting in a decay of the total population at an exponential rate. In contrast, if a re- 
flecting boundary for the Feynman integral is removed at time t = 0, the flux across the 
boundary increases as 0(\/i). If Feynman trajectories that propagate into a boundary 
are instantaneously absorbed there, the flux at such a boundary is proportional to the 
square of the local gradient of the wave function B. 

Acknowledgment: The authors wish to thank Y. Aharonov and B. Reznik for useful 
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